Modulational instability of ion-acoustic wave packets in quantum pair-ion plasmas 
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Amplitude modulation of quantum ion-acoustic waves (QIAWs) in a quantum electron-pair-ion 
plasma is studied. It is shown that the quantum coupling parameter H (being the ratio of the 
plasmonic energy density to the Fermi energy) is ultimate responsible for the modulational stability 
of QIAW packets, without which the wave becomes modulational unstable. New regimes for the 
modulational stability (MS) and instability (MI) are obtained in terms of H and the positive to 
negative ion density ratio /3. The growth rate of MI is obtained, the maximum value of which 
increases with /3 and decreases with H . The results could be important for understanding the origin 
of modulated QIAW packets in the environments of dense astrophysical objects, laboratory negative 
ion plasmas as well as for the next generation laser solid density plasma experiments. 
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In the recent years, there has been a growing interest 
in investigating various collective modes and their prop- 
erties in pair-ion plasmas (see, e.g., Refs. [T|-01 )■ Such 
plasmas are believed to be ubiquitous in most space and 
laboratory plasmas [1, Q. Moreover, it has been investi- 
gated that the pair-ion plasmas have potential applica- 
tions in the atmosphere of D-region of Earth's ionosphere, 
Earth's mesosphere, the solar atmosphere, as well as in 
the microelectronic plasma processing reactors Q- Re- 
cent investigations indicate that such pair-ion plasmas 
could also be important with regard to the diagnostic 
point of view, since the dispersion properties of wave 
modes can be used to deduce the plasma parameters 
4]. On the other hand, in view of wide applications 
in dense astrophysical environments as well as in intense 
laser produced plasmas, Misra studied the formation of 
ion-acoustic shock-like oscillations in quantum pair-ion 
plasmas |3j. Again, the propagation of wave packets in 
a dispersive nonlinear plasma medium (where the dis- 
persion is due to quantum tunneling associated with the 
Bohm potential) has been known to be subjected to the 
amplitude modulation, i.e., a slow variation of the wave 
packet's envelope due to nonlinearity. The system's evo- 
lution is then governed through the modulational insta- 
bility (MI). A number of works can be found in the litera- 
ture to study the MI in various quantum plasma systems 
(e.g., see Refs. [UTS]). 

In this work, we study the important effect of quantum 
tunneling associated with the Bohm potential on the MI 
of quantum ion-acoustic waves (QIAWs) in an electron- 
pair-ion plasma. We show that the quantum coupling 
parameter H is ultimate responsible for the stability of 
modulated QIAW packets, without which the wave be- 
comes modulational unstable. New regimes for the MI 
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are obtained with the variation of H along with the pos- 
itive to negative ion density ratio (3. We also obtain the 
growth rate of MI in terms of the system parameters. 

In what follows, we consider the nonlinear propagation 
of QIAWs in an ummagnetized quantum plasma com- 
posed of electrons and both positive and negative ions. 
The basic normalized equations read [1| 
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where the suffix a = +, — indicate the quantities for pos- 
itive and negative ions respectively. Also, <f> is the elec- 
trostatic potential normalized by ksTpe/e with ks de- 
noting the Boltzmann constant, e the elementary charge 
and Tp e = fi. 2 (37r 2 n e o) 2 / 3 /2igm e the electron Fermi tem- 
perature. Here H is the Planck's constant divided by 
2tt, m e is the electron mass and n e Q is its equilibrium 
number density. Moreover, n a denotes the ev-particle 
perturbed number density normalized by its equilibrium 
value n a o, H = hujp e /kBTF e is the nondimcnsional quan- 
tum parameter describing the ratio of the plasmonic en- 
ergy density to the Fermi energy, ui pa = \J n a Qe 2 / £ara a 
is the cv-particle plasma frequency and m a is the mass. 
The speed of the a-species particle u a is normalized by 
the quantum ion-acoustic speed c s = yTcjjX^/mT, the 
space (x) and time (t) variables are normalized by c s /oj p _ 
and uj~} respectively. In Eq. ([3]), ?(+,-) = (m,— 1) with 
to = Z + to_/Z_to + . The parameter j3 = Z+ti+qI 'Z-U—q 
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appearing in Eq. (j4]) is the positive to negative ion den- 
sity ratio with Z+ - denoting the positive (negative) ion 
charge states. The second term in the right-hand side of 
Eq. ((T|) is due to the three-dimensional (3D ) Fermi-Dirac 
pressure law for electrons given by [IJ V 
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Since the equilibrium distribution is always 3D even in 
one-dimensional (ID) geometry (we can project the 3D 
Fermi-Dirac distribution over ^-direction), the equilib- 
rium pressure must indeed be given by its 3D expression 
[Eq. ©] 0. 

In order to obtain an evolution equation describing the 
propagation of modulated QIAW envelopes we employ 
the standard reductive perturbation technique (RPT) in 
which the independent variables x and t are stretched as 
(see, e.g., Ref. [8]) £ = e(x — v g t) and r = e 2 i, where 
v g (normalized by c s ) is the wave's group velocity to be 
determined by the compatibility condition. The depen- 
dent variables (where the perturbed parts depend on the 
fast scales via the phase kx — ujt, and the slow scales only 
enter the Z-th harmonic amplitude) are expanded as 



(n a ,u o ,0 = (l,O,O) + fy £] 

n—l 1— — QO 

x exp \i{kx — ujt)l] , (6) 
where ui, k are the normalized wave frequency and wave 

(n) 

number respectively. The state variables n , , etc., sat- 



isfy the reality condition A 



(n) 



with asterisk de- 



noting the complex conjugate. Now, substituting the ex- 
pansion [Eq.©] into equations CQ)-(jll) and expressing the 
variables x and t in terms of the stretched coordinates £, 
t, and then collecting the terms in different power of e 
we obtain for n = 1, 1 = 1 the first order quantities, 
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and the following dispersion relation 
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where A = 1/3 4- H 2 k 2 /A{f3 — 1). The dispersion relation 
modified by the density ratio /3 as well as the quantum 
parameter H gives two real eigen modes for the carrier 
waves. In particular, by disregarding the negative ion 
dynamics and considering one-dimensional Fermi pres- 
sure law one can recover the previous result (l3| . In the 
short- wavelength limit (or for large k), the frequency of 
QIAWs approaches a constant value y/1 + /3m, whereas 
for large wavelength, the frequency increases with k. The 



compatibility condition is obtained from the second order 
(n = 2) reduced equations with I = 1 as 
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Proceeding in the same way as in Refs. @, [l(| and 
substituting all the above derived expressions from n = 
2,1 = 2;n = 2,1 = into the components for n = 3, 1 = 1 
of the reduced equations we obtain the following Nonlin- 
ear Schrodinger equation (NLSE) 
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^ and the dispersive and nonlinear coeffi- 
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cients P, Q are respectively given by 
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The coefficients A, B,...Ao, Bo,. ..etc., in Eqs. ([TTj) and 
(fT2"|) are given in the Appendix A. Equation (|10l) de- 
scribes the slow modulation of the first order plasma po- 
tential of QIAWs in an unmagnctized clcctron-pair-ion 
plasma. Neglecting the electron dynamics, one recov- 
ers ion-acoustic envelope solitons in a pair-ion plasma 
flij . The dispersion coefficient P = d 2 ui/2dk 2 arising 
due to quantum diffraction and the charge separation of 
the plasma particles, and the nonlinear coefficient Q, due 
to the carrier wave self-interaction, are significantly mod- 
ified by the quantum effects, as well as by the presence 
of negative ions in our plasma system. 

To study the MI of QIAW packets, we assume a 
monochromatic solution of Eq. (fT0|) to be of the form 
<fr = (f)aexp(iQ\(/)a\ 2 T) (see e.g., Refs. [1, H3]), where 
<fro is the constant amplitude of the carrier wave and 



3 



A(r) = — <5|0o| 2 is the nonlinear frequency shift. The 
amplitude and the phase of this solution is modulated 
against the linear perturbations as <f> — [4>o + 0icos(if £ — 
f2r)]exp [«Q|0o| 2r + iOicos(K^ — fir)] , where K and O 
are, respectively, the wave number and the wave fre- 
quency of modulation. We then obtain from Eq. (fpQ|) 
the dispersion relation 

n = ^Vi-§, (17) 

where K c = a/2Q/P|$o| is the critical value of K, such 
that the MI sets in for K < K c , or for wavelengths above 
the threshold A c = 2n/K c . The instability growth rate 
(letting ft — iT) is obtained as 

T = PK2 \I§- 1 > (I 8 ) 

Clearly, the maximum growth rate obtained at 
K = K C /V2 is T max = |Q||$| 2 . 

It is evident from Eq. (fTT)) that the instability condi- 
tion depends only on the sign of the product PQ. Also, 
relying on the positive or negative sign of PQ one can 
predict that when PQ > the monochromatic waves be- 
come modulationally unstable leading to the formation of 
a bright soliton, i.e. a localized pulse-like envelope mod- 
ulating the carrier wave. When P/Q < 0, the solution is 
modulationally stable, which may result to the formation 
of a dark or grey soliton, representing a localized region 
of decreased amplitude. Thus, it is reasonable to study 
the behaviors of PQ , which can be done numerically 
with the variation of the system parameters, namely, the 
density ratio /3 and the quantum coupling parameter H. 
This is the main purpose of our present work. 

We have considered the density ratio f3 to be greater 
than unity and the range of values of H < 1 for which 
the electron Fermi speed vf = \/2kBTF e /m e << c, the 
speed of light in vacuum and the quantum collective, 
mean-held effects become important (where the quan- 
tum coupling parameter gQ ~ (n e o \ 3 F )- 2 / 3 < 1, where 
Xf describes the scale length for electrostatic screening). 
The variations of P/Q with respect to k for different val- 
ues of H (with fixed m and (3) and for different j3 (with 
fixed m and H) are shown graphically in Figs. 1 and 2 
respectively. When (3 — 4, i.e., when the electron concen- 
tration is 3/4 times that of the positive ion, the stability 
(PQ < 0) and instability (PQ > 0) regions (see Fig. 1) 
for H = 0.75 and H = 0.7 are respectively < k < 0.726; 
k > 0.726 and < k < 0.4982; k > 0.4982. We find 
from Fig. 1 that for H < 0.66, PQ becomes positive 
for all values of k and any f3. The stability region de- 
creases with decreasing the values of H . Thus, decreas- 
ing the values of H or entering into the density region 
?i e o ~ 1.3490 x 10 31 m~ 3 or higher, one can find unstable 
modulated QIAWs. 

We note that in order to observe different stable and 
unstable regions for higher values of H > 0.75 we must 
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FIG. 1: (Color online) P/Q is plotted against k to show the 
stable (P/Q < 0) and unstable (P/Q > 0) regions for /3 = 4, 
m = 4 and for different values of H: H — 0.66 (solid line), 
H = 0.7 (dashed line) and H = 0.75 (dotted line). The details 
given in the text. 



have to increase the /3-values as well to satisfy the re- 
striction for quantum coupling parameter gQ as described 
above. Upon increasing the values of the density ratio (3 
keeping H = 0.7 fixed, one can find significant change of 
the stable as well as unstable regions. These can be seen 
from Fig. 2. In this figure, for f3 = 5, PQ is negative 
in < k < 0.5755 and positive otherwise, whereas for 
ft = 6, PQ < in < k < 0.64 and PQ > otherwise. 
Thus, fixing the electron concentration, and decreasing 
the positive ion concentration we can recover the sta- 
ble regions in a larger domain of k. We find that both 
the dark or grey (PQ < 0) and the bright (PQ > 0) 
soliton excitations are possible in our quantum plasma 
system. In the long- wavelength run (k < 1), the former 
dominates over a long range of values of k (for a fixed H 
and increasing while the latter exists in a bounded 
short-range of k. 

It is to be added that if in any case the plasma system 
contains negatively charged dust grains as impurities in 
the background, one can simply replace /3 by 13 — 8, where 
5 represents the percentage of negatively charged dust 
grains with respect to the negative ion concentration. In 
presence of negative dusts we may observe that for a fixed 
H = 0.7, m = 4 = /3 : PQ<0m0<k< 0.5505 and pos- 
itive otherwise for S = 0.2, and for 5 = 0.4, PQ < 
in < k < 0.464 and positive otherwise. This is how 
the charged dust impurity can modify the MI domains. 
Changing m does not modify the stability or instabil- 
ity regions, but the magnitudes of P/Q, which can be 
important for calculating soliton widths at different pa- 
rameters. The growth rate of MI is also calculated. It is 
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APPENDIX A 

The coefficients appearing in Eqs. (fTTj) and (fT2|) are 
given as follows 
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FIG. 2: (Color online) The same as in Fig. 1, but for H = 0.7, 
m = 4 and for different values of f3: /3 = 4 (solid line), /3 = 5 



(dashed line) and /3 = 6 (dotted line), 
the text. 



The details given in 



F = 



{l/3 + H 2 k 2 /fi) 
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1_ 3H 2 k 2 \ (\ H 2 ^ 
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seen that the absolute value of the maximum growth rate 
as well as the critical wave number increases (decreases) 
with increasing the /3 (H) values. 

To summarize, we have investigated the instability cri- 
teria for the amplitude modulation of QIAWs in a quan- 
tum pair-ion plasma. The quantum coupling parameter 
H is shown to play crucial roles in stabilizing the QIAW 
packets. New regimes for the MI with the system param- 
eters H and {3 are obtained in the domain of carrier wave 
numbers. The results could be important for negative ion 
plasmas, forthcoming laser produced plasmas in the lab- 
oratory as well as in dense astrophysical environments. 
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